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We present a formalism determining spin-polarized current and electrochemical potential inside 
arbitrary electric circuit within diffusive regime for parallel/antiparallel magnetic states. When 
arbitrary nano-structure is expressed by 3-dimensional (3D) electric circuit, we can determine 3D 
spin-polarized current and electrochemical potential distributions inside it. We apply this technique 
to (Cu/Co) pillar structures, where pillar is terminated either by infinitely large Cu layer, or by Cu 
wire with identical cross-sectional area as pillar itself. We found that infinitely large Cu layers work 
as a strong spin-scatterers, increasing magnitude of spin-polarized current inside the pillar twice 
and reducing spin accumulation nearly to zero. As most experimentally studied pillar structures 
are terminated by such a infinitely large layers, we propose modification of standard Valet-Fert 
formalism to simply include influence of such infinitely large layers. 

PACS numbers: 75.75.+a, 85.70.Kh, 85.70.Ay 



I. INTRODUCTION 

Spin injection, transport, and detection are key factors 
in the field of magnetoelectronics. Especially, magneti- 
zation reversal using spin-polarized current is of great 
interest 0, 0, H 0, IE1E due to its potential technologi- 
cal applications such as MRAM § , spin transistor [j| or 
spin battery [Toj . 

To understand and optimize spin-transport behavior 
in such devices, it is important to know a current dis- 
tribution in it. Particularly for MRAM applications, we 
have to know spin-current magnitude and distribution to 
optimize the current density necessary for spin-injection 
induced magnetization reversal. Variety of formalisms 
calculating magnetoelectronic transport in one dimension 
(ID) even for non-collincar magnetization has been pro- 
posed [HHm El III 

However, up to now, the spatial (3D) calculation of 
the spin-polarized current have been missing. To ob- 
tain spatial distribution of spin-polarized current (and 
spin accumulation) in a given structure, we express such 
a structure as a 3D circuit of spin-dependent-resistor- 
elements (SDRE), wherein the propagation is regarded 
as ID problem [Xulltj]. 

Resistor circuit network has been already used in ^3 
to simulate current lines in metal/insulator multilayers. 
However, in this case, they did not use spin-polarized 
current. Ichimura et al. [18| have calculated 2D distri- 
bution of spin polarized current for Co/Al lateral spin 
valve structure. They did not use the electric circuit, 
but directly solved Poison equation by means of finite 
element method. 

This article is organized as follows: Sections [H] and 
III Al present a matrix approach to calculate ID diffusive 
propagation of spin-polarized current and potential in 
the single SDRE. Section Til Bl shows how this formalism 
can be applied to a simple multilayer structure. This 



approach is just compact matrix rewriting the ID Valet- 
Fert (VF) formalism ^3 • Section IIIII explains how to 
solve general electric circuit and Section UlI Al tells how 
to divide nanostructure into circuit of SDREs. Finally, 
in Section llVl we apply our calculations to (Cu/Co) 2 , and 
(Cu/Co) 3 pillar structures, where cross-sectional areas of 
the first and last layers are assumed either the same as 
the pillar, or infinitely large. We show how the pres- 
ence of such infinitely large layers influence currents and 
spin accumulation profiles. Finally Section shows how 
to modify VF formalism to describe influence of the in- 
finitely large layers. 



II. DIFFUSIVE TRANSPORT REGIME 

In the diffusive transport regime, equations deriving 
the spatial distribution of electrochemical potential 
and spin-polarized current density j' T / inside ferro or 
non-magnetic metals are [n|, [l6|, [T^, |2(| 

V 2 (MT-Mi) = ^^, (1) 
3yi = ~^/l: (2) 

where A is spin-flip diffusion length, ayi = a(l ± /3)/2 
conductivities for up and down channels, respectively and 
e the electron charge assumed to be e = 1 in this article. 

Our model is based on the circuit of SDRE, consisting 
of layers and interfaces (Fig.^). We first express response 
of a single layer. When Eqs. are solved in ID, the 

profiles of hi and fi^ (hereafter denoted as Hyi) inside 
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FIG. 1: A general form of (a) close-end SDRE (b) open- 
end SDRE (c) classical multilayer structure having both ends 
open-ended. The arrow shows positive both structure and 
current direction. In all cases, SDRE has three layers, M — 3. 



a given material along x-axis is given by [TUllfijl 

(J*yi{x + L) = ft(x) + ^—L ± c(x) exp[— L/A] 

° a \l\ 

±d(x)— exp[L/A], (3) 
a \li 

where j c h = j\ +j\ is charge current density and positive 
current direction is towards positive rr-direction. Ener- 
gies c(x) and d(x) are amplitudes of exponential dump- 
ing of /if /I and the energy ft{x) is the asymptotic elec- 
trochemical potential equivalent to weighted average of 
electrochemical potentials ft = {<jyi-\ +a] r ^i] r )/u. The en- 
ergies ft, c and d will be determined later from boundary 
conditions. 

The relation between [tyy jyi and ft, c, d used in 
Eqs. J2| and © can be expressed in a compact matrix 
form 
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where J 
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sj 



yi is up/down spin-polarized current 
with cross-section area S of a given layer, J c h = Jf + Jj 
is a charge current flowing through layer. The scaling 
factor it which has no physical meaning, is introduced 
to adjust units of newly defined H and F vectors. For 
numerical calculations, value of u should adjust order of 
fiy I and uJ-^n. When expressing any measurable quan- 



tities (e.g. /i T/i , 
depend on u. 



Jyi, resitances), they obviously do not 



The left side of Eq. (J3J (.ff -vector) contains variables, 
which are conserved across interfaces or electric nodes. 
On the other hand, the right side of Eq. Q (.F-vector) 
contains variables, which are used to calculate propaga- 
tion of the Jyi and Hyi through layers [Eq. ©]. Hence, 
the dynamic D-matrix relates coefficients at the bound- 
ary to coefficients for propagation. This approach is well- 
known in the optics of anisotropic media |2lL . 

Propagation of .F-vector through layer with the thick- 
ness ( "length" ) L is expressed by a propagation P-matrix 
[Eq.©] 
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F x+L = [P]--F x . 

Hence, the relation between H = \py fj,^, uJy uJ^ 
both ends of the layer is expressed as 

H x = DPD 1 H x+ l = KH x+ l. 



(5) 
at 

(6) 



A. Interface resistivity and shunting interface 
resistance 

In the previous Section, we have expressed response 
of one layer in homogeneous material. Here, we de- 
scribe the interfacial properties including (i) interface re- 
sistance ARyi — 2AR*(1 qp 7) and (ii) interfacial 
spin-flip scattering, described by a shunting resistance 
AR S , "shortcutting" up and down channels at the inter- 
face. The response of interface can be expressed by Kj nt 
matrix 
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Notice that alternative interfacia l sp in-flip scattering by 
oj-parameter was introduced in [23, |24( describing spin 
relaxation at the interfaces by a thin interfacial layer of 
spin- flip- length A/ = tj/5j where tj is interfacial layer 
thickness. From comparison of K-matrices from Eqs.® 
and ©, AR S = 4t] /(a 16 j)sinh Si where tj and 07 are 
interfacial layer thickness and conductivity, respectively. 
When interfacial spin-flip scattering does not exist, then 
i?, = 00. 



B. Simple multilayer structure 

Although the formalism developed here is mainly to 
calculate jJtyy jyi inside an electrical circuit, we shall 
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first show how presented 4x4 matrix algebra can be used 
to calculate electrical response of a multilayer structure 
[Fig.Hc)]. _ 

The description for a single layer by matrix K in Eq. © 
gives relationship between fj^ n and Jf n at both ends 
of the layer. Hereafter, instead of continuity of spin- 
polarized current density j| / j over interface, we consider 
continuity of spin-polarized current Jyy to take into ac- 
count variable cross-sectional area S of layers. The rea- 
sons and validity are discussed later in Section^ Ob- 
viously, when all layers has the same S, our formalism 
provides the same results as original VF formalism [llj . 

The boundary conditions at interfaces are continuity of 
/if/I and Jyi, i.e. continuity of if- vector. Consequently, 
the response of whole multilayer structure can be written 
as [21I, 



p(0) = M (M+1) , F (M+1) 



(9) 



where 

m( m+1 > 



1 K&> K« KB 



K (M-1) K (M) K (M) D (M+1) (1Q) 



where upper index in parenthesis denotes for interface or 
layer number, M the number of layers. 

Because materials (0) and (M + 1) are semi-infinite, 
/if n inside them must not exponentially increase. Hence, 
some exponential terms in Eq. must vanish, namely 
c<°) = 0, d< M+1 ) = 0. Hence, the vectors F (0) and 
p(M+i) are ii m ited f orm 
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Now, Eq.© can be solved. Substituting Eq. l(TT|) to 
Eq. (jHJl, we found all unknowns in i^( M+1 ) 
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(12) 

where My are elements of M-matrix. The value of /x^ - 1 
can be arbitrary value, as it only adds a constant to the 
profiles of fXyy Other simplification of Eq. I|12|) follow 
from M11 = 1 and M21 = 0, in consequence of (i) for 
Jch — the terms fJ-yy must be constant and equal 
to each other and (ii) for J ch = 0, C ( M+1 ) = 0. Then, 
solution of Eq. (|12[) can be written as 
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Now, /2( M+1 ) and c' M+1 ' are known and therefore the 
vector F ( M+1 ) can be reconstructed from Eq. i|ll|) . Con- 
sequently, the profiles of ^yy fi, Jyi etc. in the entire 



structure can be determined by recursive applying step- 
by-step matrix multiplications in Eq. (|1(J|) 

Finally, resistivity of the multilayer structure (between 
first and last interface) is R = {ft M+v > - ^°y/J ch . 



III. ELECTRIC CIRCUIT 

As demonstrated in the previous Section, the M- 
matrix [defined by Eq. can characterize the entire 
multilayer structure where the same charge current J c h 
flows across all layers. In this Section, we extend the 
previous formalism to an electrical circuit (network) of 
SDREs, mutually connected at nodes. In general, the 
SDRE is composed of any sequence of materials (layers) 
and/or interfaces, as depicted on Fig. ^ There are three 
types of SDRE, depending whether the length of SDRE 
is finite or infinite: 

• close-end SDRE [Fig. ^a)] has finite length and 
hence both ends of this SDRE are attached to 
nodes. Because the boundary condition at nodes 
are described by if-vectors, whole SDRE is de- 
scribed by K' b ' -matrix relating H- vectors at both 
ends of SDRE: 
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where, for later purpose, can be rewritten into 
four 2x2 submatrices 



(15) 



where 6 denotes a SDRE number in the circuit 
and M;, is the number of layers in 6-th SDRE, 

tlMM b ) = [^K^^KM^T and j[b](M b ) = 

^j[b]{M b ) ^ j[b](M b )jj: ^ t denoting vector transposi- 
tion. Hereafter, indices in square bracket de- 
notes for SDRE number, whereas indices in ordi- 
nary parentheses denotes for number of layer inside 
SDRE. Positive current direction is from layer (1) 
to layer (M b ). Analogous to Eq. (jTU|) . consists 
of layer and interface contributions 



[b] _ K [b](i) kWW 



K 



hit 



K 



[b](M b -l) K [b](M b ) 



(16) 



When SDRE contains only one layer (i.e. it consists 
of single metal), then [Eq. ®] K [fc] = K w(1) = 

D [6]p[61[ D [61]-l_ 

open-end SDRE [Fig.^b)] has one end "finite" and 
connected to the node. The other end is infinitely 
long and at its end the charge current J^°\ flow- 
ing into the 6-th open-end SDRE, is applied. Be- 
cause boundary conditions on node are described 
by iJ-vector and boundary conditions of the prop- 
agation into infinity by _F-vector, the 6-th open-end 



SDRE is described by Z' '-matrix as 

pMO) _ z [b] H [fe](M i) ) j / 17 n 

where we have used convention that direction of 
SDRE and positive current direction goes from in- 
finite end of SDRE toward its finite end [Fig.[IJb)]. 
Analogous to Eqs. fTHjl and lfTfi)l 

Z [b] = [ D M(0)]-1 K |6](0) K [6](1) K [6](l) 

K M(M 6 -l) K [6](M i ,)_ ^ 

When the SDRE contains no layer (i.e. it contains 
only single material continuous to infinity), then 

Z [b] = [D [6]]-1, 

• multilayer structure [Fig. uJL c )] described in 
Sect. Ill Bl may be understood as a special type of 
SDRE, having both ends open-ended. 

In general, the electric circuit is assumed to have N 
nodes, C close-end SDRE and E open-end SDRE. 

The boundary conditions valid for each node are de- 
termined by generalized Kirchoff's laws: 

mL 6 "t/4 = const «,T/l = Mn.T/ii ( 19 ) 

E J M/1 = °. (2°) 

i.e. the values of /z-j- /^ has to be identical for all ends of 
SDRE connected to each node and sum of polarized cur- 
rents J| /i entering each node has to be zero. Subscript n 
denotes for node number, n = 1 . . . N and b n = 1, . . . B n 
denotes the SDREs connected to the n-th node. 

Hereafter, we use two u-notation (i) £t ni -m relating 

/if/I at the n-th node and (ii) ^yf\ f^y^ b ^ denoting 
fly I at start, end of the 6-th SDRE, respectively. The 



relation between these two /z-notations is given by con- 
nections between nodes and SDREs, i.e. by the topology 
of the electric circuit . 

Following the previous discussion, the problem is how 
to treat large number of linear equations lfTI)l JT7|| l (T5)l (f^Of) 
giving relations between /Lt-p n and Jj n on the nodes and 
between ends of SDREs. We do it by solving one large 
matrix expression, 

Q • H = I, (21) 

to which we apply all the previously mentioned rules. 
The H-vector contains (J> n ,t/l f° r all nodes and also 
uJ^J^, b = 1 . . .C + E, at the end (i.e. at the M b - 
side) of each SDRE. Hence, H-vcctor has 2(C + E + N) 
elements. 

The role of the Q-matrix consists of three parts 
[Eq.|23]: 

(1) to relate [tyi, Jy± at the "finite" end of the open- 
end SDREs and between boundary conditions of 
propagation towards infinity. This is provided by 
2E linear equation, and hence this part occupies 
2E rows in Q-matrix. 

(2) to relate ny^ between start and end of the close- 
end SDREs (2C linear equations). 

(3) to realize current conservation at each node (2N 
linear equations). 

These contributions are studied in detail in following. In 
total, Q-matrix has 2(C+E+N) rows, the same as length 
of JJ-vector. So Q is a square matrix. 

An example of a circuit of SDREs is depicted in Fig- 
ure El The circuit consist of N = 3 nodes, connected 
by SDRE. The circuit consists of C = 3 close-end SDRE 
and E = 2 open-end SDREs. Then, the resulting equa- 
tion QH = F looks like: 
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(22) 



where fi n = [jU n ,T > f i n,i] T and 1, denote 2x2 unitary matrix in Eq. IL'21) represents two rows (for up and down 
zero matrix, respectively. Each row which is shown in Q- channel), so hereafter we call it "double-row". 
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FIG. 2: The schema of a circuit example, consisting of three node, two open-end SDRE and three close-end SDRE. The arrows 
parallel to SDRE shows the defined direction of SDRE (and hence also positive current direction). 



(ad 1) The role of the first part in the Q-matrix ond double- row in Eq 
relates both ends of open-end SDRE, described by 
Eq. CO, = Z^H^ Mb \ However, we 

only know 2 variables out of 4 in vector F ^'°- ) = 
^M(o) iC [6](o) )d [6](o) )Ui/ W(o)]T, We know charge current p'[b'K0)_ 

J^ I(0) entering 6-th SDRE and that c [bI(0) = 0, because 
/i| n must not exponentially increase towards infinity 

[Eq. ifTTl) ]. Hence, taking from Z' b '-matrix [Ea. l(T7|) ] only 
rows corresponding with known value in F^- b ^°\ we get 
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This equation is substituted into first part of in Q-matrix 
[Eq. (J23)] in the form 



z'V 1(Mi> 



fb](o) 

In the case of the last open-end SDRE connected to 
the structure (in our example SDRE number b' = 5), the 
situation is slightly different. In the vector f[ 6 ]( ) = 

[^'](0) )C [6'](0) )d [b'](0) !U j[6'](0)p (a) we have tQ get a 

value of fi\ b ^°\ which sets an absolute value of all fi^n 
and y, inside circuit to an arbitrary value and (b) getting 

a charge current J^}^ would be redundant as a sum of 
all charge currents entering circuit has to be zero. Hence, 
the last double-row in the first part of the Q-matrix (sec- 



and is substituted to Q-matrix analogous as Eq. 

(ad 2) The second part in the Q-matrix gives rela- 
at the start of the 6-th close-end 



jfj{ Mh) at the end of 6-th SDRE. 



tion between fiy^ 

SDRE and ^}f h \ 
This relation is given by the first double row taken 
from Eq. (|15p. It is substituted into Q in the form 

,M(0) + K [b] [b](M b ) + uK [b] 7 J [hl(Mi>) . 



= 

(ad 3) The last part in Q-matrix describes the cur- 
rent conservation at each node, described by Kirchoff's 
law Eq. (|20|) . As specified, H- vector contains only values 



of current at the end of each SDRE, J. 



[6] (Mi, 



Hence, if 



the 6-th SDRE starts at the n-th node (as No. 2 and 4 
SDREs at the node 1), then the current is expressed 
by the second double-row in Eq. from / u f/i Mi ' ) ' 

j\ b )l Mb) as ujWW = K ™/iMW +uK [ fjJ^ Mb \ For 
example, in case of node n = 1, the current con- 
servation is -jf}™ + jf}® + jf}® = 0, which 
is substituted to Q-matrix to the 6-th double-row as 



0. 
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FIG. 4: j sp through (Cu/Co) 2 structure with constant cross- 
sectional area S calculated by ID model (line), by our 3D 
model with (x) and without (o) A-renormalization. 



FIG. 3: An example of dividing nanostructure into 3D circuit 
of SDRE. Each wire represents a "bus" of channel-up and 
channel-down. The inset remark that each SDRE consists 
of resistivity for channel-up, for channel-down and a spin- 
flip-scattering resistance between both channels. The large 
resistors denotes for interface resistivity [Eq. JHJ]. 



Although the construction of the Q-matrix as pre- 
sented here may be tedious, it is rather direct to establish 
its construction numerically. When the equation QH = F 



is solved, values of fJiyi, J-m for each SDRE are directly 
written in H- vector; their profiles can be determined by 
step- by step applying multiplication in Eas. (|14()(|16() and 

©Gil!. 



A. Construction of 3D electric circuit 

In the previous Section we have derived the formalism 
to calculate J-f n and fij n in arbitrary electric circuit. In 
this Section, we explain, how to describe the whole nano- 
structure as a circuit of SDREs, as sketched in Fig. 01 

Each part of the nanostructure is divided into 3D rect- 
angular grid, the circuit nodes positions being Xi, yj, 
Zfe. Then, a given SDRE, for example in the x-direction, 
has length L, t = (xi-i + Xi)/2 and cross-sectional area 
Sjk = UjZk- Due to this treatment, the grid does not 
need to be equally spaced, but just rectangular. 

At the interface between two different materials, 
e.g. A and B (see Fig. |3J), SDRE is described by 
Ka^b matrix consisting of three contributions Ka^b = 
Ky^intKintKint^s [Eq. (TgJ]. The K^int is contribu- 
tion from grid point (node) inside A material to the in- 
terface with (SDRE is now along the z-direction) — 
z-mt~ Zk-i and Sij — Xi]jj. K lnt is the interface resistivity 
matrix given by Eq. (JSJ . Kj nt ^B is contribution from the 
interface to the node in B material having L — — Zmt 
and Sij = Xiyj. 



SDRE circuit the volume of metal in all resistors is three 
(two) times larger than in reality. In such a case, the 
spin-polarized current would have larger dumping, as it 
would diffuse into larger volume. To correct this, it is 
necessary to calculate ID propagation of /if n in each 
SDRE by slightly modified Eq. JTJ 



/ 



<9 2 (MT ~ Ml) _ MT - Ml 



dx 2 



X 2 



(25) 



The electrical circuit network as described above de- 
scribes correctly charge current. However, inside 3D (2D) 



where / is a dimension of the SDRE circuit. In other 
words, when a given nanostructure is described by 3D 
(2D) electrical circuit, the A should be increased by a 
factor V3 (\/2). This A-normalization should not be ap- 
plied for open-end SDREs, as their contribution is cor- 
rectly described by ID propagation. 

The advantage of this A-normalization is shown in Fig- 
ure 0J where j sp through (Cu/Co) 2 structure (described 
and studied in detail in the next Section) is compared 
between ID model (full line) , and our 3D model without 
A-normalization (o) and with A-normalization (x). We 
can see that with (without) A-normalization, the agree- 
ment between ID and 3D model is about 1% (30%). The 
same is valid for spin-accumulation Afi = /x-j- — /i^, where 
calculation without A-normalization leads to 30% smaller 



B. Surface scattering 

In this Section, we describe how to incorporate surface 
scattering to the presented formalism. 

Surface scattering can be described by a shunting resis- 
tance R s shortcutting up and down channels for a nodes 
situated just close to the wire surface. R s has value 
R Si n — AR S /S n , AR S being surface scattering resistiv- 
ity and S n being surface area corresponding to the n-th 
node. When surface scattering is not presented, then ob- 
viously R s .n — inf. Surface scattering can be described 
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FIG. 5: Sketches of studied type structures, (a) S-constant (infinitely long nanowire with constant cross-section area S), (b) 
column (infinitely long pillar deposited on infinitely large Cul layers) and (c) constriction (both Cu3 cover and Cul buffer 
layers are infinitely large) (d) constriction with Col infinite layer. 
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FIG. 6: j sp along center axis of the (a) (Co/Cu) 2 TT, Tl (b) (Co/Cu) 3 TTT, TIT (c) (Co/Cu) 3 TTI, III for S-constant, column 
and constriction structure types [Fig. [5]. 



by modification of Kirchoff law [Eq. l|20jl] 

[k n ] _,_ Mn,T ~ 



2^ J n,m + 



Rs 



0. 



(26) 



To incorporate this modification into Q-matrix for the 
?i-th node, we add G s matrix 



-1 
1 



(27) 



on the n-th double-column and ro-th double-row in the 
last part of Q-matrix which describes current conserva- 
tion in each node. In our example given by Eq. ill' 21) and 
Figure |3 to add R s to n = 2 node, we place G s to the 
7-th double-row and 2-nd double-column into Q-matrix. 



IV. APPLICATION TO (Cu/Co) PILLAR 
STRUCTURE 

In this Section we use the above developed formalism 
on (Cu/Co) 2 and (Cu/Co) 3 pillar structures. We will 
show how spin-polarized current density j sp = — jy 
and electrochemical potential ^ vary differently be- 
tween when whole structure is a pillar or only a part of 
the structure is a pillar attached to an infinitely large 
continuous layer. In the literature, pillars terminated 
with infinitely large continuous layers are commonly used 



We have studied structures consisting of 2 and 3 
Co layers, called (Cu/Co) 2 , and (Cu/Co) 3 with di- 
mensions in nm Cul/Col(40)/Cu2(6)/Co2(2)/Cu3 and 
Cul/Col (40) /Cu2(6) /Co2(2) /Cu3(6) /Co3(20) /Cu4, re- 
spectively. The square-shaped pillar 100 nm in size be- 
gins with the Col/Cu2 interface. The considered struc- 
ture types are defined in Fig. [S] as (a) cross-sectional 
area S constant, (b) column, (c) constriction and (d) 
constriction with infinite Col layer. In the following 
discussion, "constriction" corresponds to the case (c). 
The "infinite" homogeneous layers were approximated 
as a square pillar of 800 nm in size. The magnetization 
of Col layer is always fixed as "up" (|), whereas mag- 
netic orientations of other Co layers are varied. Note 
that in our diffusive transport calculations, the magneti- 
zation orientation with respect to the structure (e.g. in- 
plane or out-of-plane) do not play any role, but only mu- 
tual magnetization orientation (parallel or antiparallcl) 
does play an important role. The charge current passing 
through structure is assumed to be J c h = 1 mA, equiva- 
lent to the averaged charge current density in the pillar 
jch = 100xl0 9 A/m 2 . The electrical properties of materi- 
als are of room temperature 0, 153, : electric conduc- 
tivity ctcu = 4.81xl0 6 fi- 1 m- 1 , aco = 4.2xl0 6 fr 1 m- 1 , 
spin-flip-lcngths Aq u = 350 nm, Aco = 60 nm and Co spin 
bulk assymetry 7 = 0.35. We assume no interface resis- 
tance and no interface and surface scattering. The SDRE 
grid size is 10 nm. 
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j Bp [10 9 A/m 2 ] 




meV] 


MR [%] 


(Cu/Co 


2 It S-constant (+) 


14.85 [14.93] 


-0.223 


[-0.220] 


0.483 [0.485] 


(Cu/Co 


2 It column (□) 


19.0 [19.8] 


-0.267 


[-0.276] 


0.63 [0.66] 


(Cu/Co 


2 || constriction (O) 


29.4 [33.7] 


-0.052 


[-0.009] 


1.01 [1.17] 


(Cu/Co 


2 || constr, Col infinite 


28.7 [23.3] 


-0.050 


[0.001] 


1.43 [2.48] 


(Cu/Co 


3 HI S-constant (x) 


19.51 [19.67] 


-0.078 


[-0.077] 


0.444 [0.448] 


(Cu/Co 


3 HI column (A) 


23.6 [24.5] 


-0.171 


[-0.184] 


0.55 [0.57] 


(T'ti /rv. 

^V_jU/ V_jO 


constriction ( ^) 




n mo 


[-U.UUOJ 


U. / [U.olJ 


(Cu/Co 


3 HI S-constant (*) 


5.94 [5.93] 


-0.442 


[-0.438] 


0.116 [0.115] 


(Cu/Co 


3 HI column (>) 


8.08 [8.49] 


-0.490 


[-0.494] 


0.162 [0.171] 


(Cu/Co 


3 TTI constriction (<) 


9.7 [10.5] 


-0.458 


[-0.456] 


0.192 [0.210] 



TABLE I: Values of j Bp , A/i = /U| — /^j at the position of free Co2 layer for (Cu/Co) 2 and (Cu/Co) 3 structures. Magnetoresistance 
ration (MR) is determined between first and last Co/Cu interface. All values are determined as averadge over whole pillar 
cross-section area. In square brackets we present values calculated from modified ID VF formalism, which takes into account 
a variable cross-sectional area of the layers (Sect.[VJ. Symbols in parentheses denote structure notation in Fig. 1111 



A. Current density in the structure 

Figurc|Slshows the profile of j sp along center axis of the 
structures: (a) (Cu/Co) 2 (b) (Cu/Co) 3 with parallel Col 
and Co3 layers and (c) (Cu/Co) 3 with antiparallel Col 
and Co3 layers for S constant, column and constriction 
type structures. In all cases, j sp for parallel Col and 
Col layer is larger than for antiparallel configuration. 
Furthermore, j sp is enhanced at the position of free Co2 
layer for the column and constriction types compared 
to the S constant type structure. For example, in the 
case of the constriction type structure with (Cu/Co) 2 || 
configuration, j sp is enhanced by a factor of 1.75 and in 
the case of (Cu/Co) 3 |H by 1.5. When only Cu buffer 
layer is infinite in size (column-type structure) , j sp at the 
position of Co2 layer is enhanced too, but not so strongly. 
Figure HJc) shows that for (Cu/Co) 3 with antiparallel 
Col and Co3 layers, j sp is significantly reduced for all 
types of structures. 

The origin of the j sp enhancement is following: the 
pillar is attached to infinitely large Cu layer which pro- 
vides large volume for spin-current to be scattered and 
so acting as a strong spin-flip-scatterer. Hence, infinitely 
large Cu layer works as a small shortcutting resistance 
between up and down channels. Consequently, short- 
cutting of up and down channels leads to an increase 
in j S p. The increase of j sp is related with increase of 
spin-polarization efficiency p = j S p(^ch/5'piiiar), as charge 
current flowing though pillar J c h is fixed in all our cal- 
culations. Consequently, increase of p leads to decrease 
of critical switching current J SjC h which is necessary to 
reverse magnetization direction of free layer. 

As the constriction type structure is in the most com- 
mon use for Co/Cu pillar structures, the j sp enhancement 
(i.e. J^ch reduction) has been already widely used 0, Q 
without being noticed. Similar effect can be realized by 
inserting a layer with small characteristic spin-flip resis- 



tance AR\ = X/a, such as Pt, Ag, Au, Ru above the last 
Co layer or bellow the first Co layer. Such cover layers 
has been used [2^, |3(| since the first pioneering experi- 
ment H , but their contributions to the j sp enhancement 
have been observed recently [3(], El 113 • 

Figure shows profiles of charge current j c h = jf + j[ 
along the center axis of the (Cu/Co) 3 structure for S 
constant, column and constriction type structures. Ob- 
viously, for S constant structure, the j c h is constant. In 
the case of infinite Cu termination, j sp decreases approx- 
imately exponentially over the characteristic length of 
50 nm. The same decay of j sp is presented in Fig. [5] for 
infinitely large Cu layers. 

In Table [0 we summarize averadged values of j sp in 
the position of Co2 layer in all the types of studied struc- 
tures. These j sp values may differ from those presented 
in Fig. due to lateral inhomogenity of j sp inside pil- 
lar, discussed in next Section. The largest averaged j sp 
is obtained for (Cu/Co) 3 |TT constriction structuture 
(31.0 x 10 9 A/m 2 ) and (Cu/Co) 2 ft constriction struc- 




i . i . i . i . i ii. i . i . i . i 



60 -40 -20 20 40 60 80 100 120 140 
distance from pillar foot [nm] 

FIG. 7: jch along center axis of the (Co/Cu) 2 structure for 
S-constant, column and constriction structure types [Fig. El- 
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FIG. 8: Spin polarized current in (Cu/Co) 2 constriction-type structure (a) without infinitely large Col layer, TT (b) with 
infinitely large Col layer TT (c) as (b) but TJ,. 
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FIG. 9: (a) j sp and (b) j c h in the middle of the Co2 layer inside (Cu/Co) 2 structures. 



ture (29.4 x 10 9 A/m 2 ), providing j sp enhancement by 
a factor 2 with respect to (Co/Cu) 2 S-constant structure 
(14.85x 10 9 A/m 2 ). This tendency is already explained 
in above paragraph. 

B. Current inhomogenity inside pillar 

Figure [S] presents a map of spin-polarized current den- 
sity j S p inside (Cu/Co) 2 constriction type structure for 
(a)(b) parallel and (c) antiparallel magnetization config- 
uration. The case (a) is a structure without infinitely 
large Col layer whereas (b)(c) with it. A map of j c h is 
not presented here as it looks similar to j sp with TT. The 
jsp inside Co2 layers is more homogenous and flows well 
perpendicular to the interfaces although in the surround- 
ing Cu layers the j sp can have rather large inclination 
and inhomogenity. The j sp tends to be more homoge- 
neous when passing Co layer, causing that the j sp in the 
adjacent Cu layers can have rather large in-plane com- 
ponents. This is remarkable in the case (c). Very similar 
tendency is found for j c h- 

The above mentioned characteristics of j sp are 
consequence of larger spin-flip resistance AR\ = 



A/er of Co with respect to Cu, Ai?ACo=14.3 film 2 , 
AR\ t cu=7 -3 fi^m 2 . It means that spin-flip is more likely 
to occur inside Cu than inside Co. In other words, Co 
is "harder" material than Cu for j sp to penetrate into 
it. The above mentioned characteristics of j c h are simply 
consequence of ceo <C <7c u - 

Figure shows cross-sectional (in ^-direction) profiles 
of (a) jsp and (b) j c h through pillar in the middle of the 
Co2 layer for (Cu/ Co) 2 structure. Both j sp and j c h are in- 
homogenous having minima at the structure center. This 
is due to inhomogcnous current injection into the pillar 
from infinitely large layers. There is no j S p and j cn in- 
homogenity for S constant and nearly no inhomogenity 
in case of column type structure as the j sp and j c h are 
homogenized by a Col layer. For constriction type struc- 
ture the inhomogenity is 12% for TT and 5% for TJ.. Inho- 
mogenity is defined as ( j max - j m i B )/ (jmax + jmin), where 
maximal, minimal value is taken from cross-sectional cur- 
rent distribution in Fig. |5J The inhomogenity is reduced 
for 11 due to different conductivities of up and down 
channel conductivities, leading to other current homog- 
enization in Cu2 spacer layer. If Col layer is infinitly 
large, the inhomogenity is increased to 10% and 29% for 
|'T and TJ,, respectively. The j cfl inhomogenity is 8% and 
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FIG. 10: Electro-chemical potential fi yi along center axis of the (a) (Cu/Co) 2 TT, Tl (b) (Cu/Co) 3 |Tt> TIT and (c) (Cu/Co) 3 
TTTi TTT f° r S constant, column and constriction structure types [Fig.Q3- 



20% for constriction type structure with and without in- 
finitely large Col layer, respectively. 

Here remains a question whether such a j sp inhomogen- 
ities make current magnetization reversal easier or not. 
Adventage of j sp inhomogeneity may be, that it localy 
enhanced j sp near the pillar edge, wheas decreasing j sp 
at the pillar center. As shown in TableQJ the mean value 
of j S p is about the same for (Cu/Co) 2 constriction struc- 
ture with and without infinitely large Col layer. How- 
ever, in the second case, j sp is much more inhomogenous. 
This is not well presented in Fig.[5Ja), we do not see the 
largest j sp flowing in the vicinity of the corners of the 
square pillars. Disadventage of j sp inhomogenity may be 
different magnitude of torgue exerted on magnetic spins 
of free Co2 layer. This may be particulary important in 
the case of high speed switching associated with magne- 
tization precession. 

C. Electrochemical potential inside structure 

Figure ITTTI presents profiles of /if, /if and fb (hereafter 
/i-profiles) along the center axis of the (a) (Cu/Co) 2 TT, 
U (b) (Cu/Co) 3 TTT, TIT and (c) (Cu/Co) 3 HI, TU- 
S-constant structure is presented for both magnetization 
directions of Co2 layers, although column and constric- 
tion types are presented only for f magnetization of Co2 
layer. In contrast to j sp , inhomogenity of spin accumu- 
lation A/i at the position of free Co2 layer is very small, 
mostly bellow 1%. Table H] presents a mean value of A/i 
for all types of the studied structures. 

Figures EH ( a )(b) and (c) show that /i-profiles de- 
pend slightly on magnetization of free Co2 layer because 
tco2 "C Aco- Figure ITUT a) exhibits suppression of A/i in 
the vicinity of infinitely large Cu layer. The reason is 
exactly the same as discussed in Sec. II V A1 the infinitely 
large Cu layer works as a strong spin-scatterer, causing 
a small spin-flip resistance (large scattering) between up 
and down channels. Obviously, such a shortcut reduces 
A/i. 

This is contradictory to [23, On], [22 , where it is argued 
that presence of spin-scatterer increases spin acululation 



A/i inside pillar. It should be emphesised that presence 
of spin-scatterers increases j sp (and magnetoresistance) 
in the pillar, but reduce A/i. 

Table [I] shows that the largest A/i is obtained for col- 
umn type structure, by 20% larger than for S-constant 
type structure. The reason is as follow: when Cul is not 
infinitely large (S constant structure), A/i changes its 
sign approximately in the middle of Col layer [Fig JlOf a) 
S constant]. When Cul is infinitely large, it acts as 
strong spin-scatterer and shortcuts up and down chan- 
nels. Hence, A/i at Cul/Col interface is nearly zero and 
hence larger A/t is obtained at the Col/Cu2 interface 
and inside Co2 layer [Fig. HOf a) column]. But to realize 
this, it is necessary that up and down channels above free 
Co2 layer should not be shortcut meaning that the cover 
layer has not to be infinite. 

Figure HDT bl shows that A/i is reduced when Col 
and Co3 layers have parallel magnetization configura- 
tion. Figure HOf c) shows an increase in A/i when Col 
and Co3 layers are antiparallel, enhancing A/i by factor 
of 2 with respect to (Cu/Co) 2 S constant structure. In 
this case the types of structure are not so important as 
in (Cu/Co) 2 case, because A/t at the position of Co2 free 
layer is "screened" by spin-scattering inside Col and Co3 
layers. 

Hence, there is effectively no spin- accumulation A/i in 
the case of a commonly used constriction type structure 
with two FM layers. It may explain why this contri- 
bution to magnetization reversal (predicted in |33() was 
not observed in Co/Cu structure |25| . To obtain non- 
zero A/i, it is necessary to use either 3 FM layer system 
with antiparallel configuration of first and last FM lay- 
ers, or to ensure that structure above free layer will not 
contain spin-scatterers. It can be reached when pillar 
structure above free layer (i) does not contain any strong 
spin-scatterer layers (as Au, Ag, Pt, Ru) and (ii) the 
cover layer does not contain large volume of metal. It 
means that pillar current drain should be realized by a 
long pillar or by a thin cover layer. 
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FIG. 11: Dependence of (a) j Bp and (b) A/i on AAR for different studied (Cu/Co) structure, symbol notation in Tab|I] Lines 
are the best linear fits, dashed line and bold symbols (full line and normal symbols) denote parallel (antiparallel) Col and Co2 
layers. Cul infinite means column and constriction types, Cul finite means S constant type. 



D. Magnetoresistance 

Finally, we discuss influence of type structure on mag- 
netoresistance ratio (MR), presented in the last column 
in Table [I] The value of MR is determined as MR= 
(A/if — A/i|)/(A/i| + A/2;), where |, j denotes magne- 
tization of free Co2 layer, Ap,y l = /*i as t,t/j. - MfhsU/i, 
where jSfirst,T/l an( i Aiast,T/l are determined on Cu side 
of the first and last Cu/Co interface, respectively. 

The calculated value of MR in the case of (Cu/Co) 2 
S-constant structure is 0.48%. However, in the case of 
(Cu/Co) 2 constriction type, it reachs 1.01% (enhance- 
ment by a factor of 2) and in the case of the constric- 
tion with infinite Col layer even 1.43% (enhancement by 
a factor of 3). The last value may be misleading as this 
increase is mainly due to resistivity reduction of Col in- 
finite layer. 

Table [i] shows MR is strongly affected by the type 
structure, increasing with an increase of j sp rather than 
A/i. This tendency may be observed also for (Co/Cu) 3 
layer, showing increase (decrease) in MR when Col and 
Co3 layers have parallel (antiparallel) magnetization con- 
figuration. 

Figure ITlT a) shows dependance of j sp on AAR, where 
AR, = i?| — where |j I means up, down magneti- 
zation of Co2 layer, respectively. It shows that j sp is 
proportional to AR. This figure is analogous to the ex- 
perimental dependence of critical switching current J s , c h 
as a function of AR |3(J, 13 ll l32l | , where they found that 
l/J s ,ch is proportional to AR. 

These graphs are analogous from following reasons: to 
switch magnetic layer we need to overcome the critical 
spin-polarized current density 

js,sp,0 = pJ s ,ch/ S p il\ al , (28) 

where J S;C h is critical switching charge current and p 
is spin-polarized current injection efficiency. Because 
js,sp,o an d Spiiiar are constant, so p ~ l/J SjC h- On the 
other hand, from our calculations, we can express p as 
P = jsp/(Jch/ S 'piiiar), where J c h is a fixed charge current 



flowing to the structure; therefore p ~ j sp . However, in 
our model, j sp is slightly larger for Col and Co2 parallel 
configuration (Fig. [HJ), although in j33, 0, H2, ^/J s ,ch 
is larger for antiparallel configuration, in agreement with 
the spin-transfer model |34|. Figure ITlT a) also provides 
an important consequence: when changing electric prop- 
erties of surroundings of FM(fixed)/spacer/FM(free) lay- 
ers, an enhancement of j sp in the position of free layer 
leads to enhancement of MR. 

Figure UlT b) shows dependence of A/i on AR. We can 
see that with increasing AR, A/t is reduced. Two dif- 
ferent slopes correspond to two different "sources" A/t 
acting as with different "hardness" . The source hard- 
ness is determined by a presence or absence of scatter- 
ers bellow fixed Col layer. Hard source is for constric- 
tion and column structure type, i.e. when Cul is in- 
finitely large. Weak one is for S constant structure type, 
i.e. when for Cul is not infinitely large. The explana- 
tion of this behavior has been already provided in pre- 
vious Section llV CI when spin- flip scatterer is presented 
bellow fixed Col layer, it vanishes A/i on Cul/Col in- 
terface and hence provides harder source of A/t. Fig- 
ure 1111 also shows that when changing surroundings of 
FM(fixed)/spacer/FM(free) layers, an increase in A/i at 
the position of free layer is related to a decrease of MR. 

Remark, that for constriction type structure with Col 
layer infinitely large [case(d) in Fig. [5], the "source" of 
A/t becomes softer than for constriction type [case (c) in 

Fig.ua. 



V. EXPRESSION OF INFINITELY LARGE 
LAYERS BY ID MODELS 

As we have shown, the constriction type structures 
modify the profile of both j sp and /tf^ with respect to S 
constant types, which provides equivalent results to ID 
VF formalism pd|. However, to describe j sp and /t-m 
inside structure by 3D model may be tedious procedure. 
That is why here we discuss briefly, how to modified VF 
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FIG. 12: (a) averaged value of spin-polarized current and (b) Hf/i in the structure center calculated for constriction type 
structure with and without Col infinite layer in the (Cu/Co) 2 Tl compared with ID calculations taking into account infinite 
Cu layers. See Section [V] for detail. 



model to take into account infinitely large layers. 

We propose, that each layer can have its own cross- 
sectional area Si. Then, instead of conserving spin po- 
larized current density j| / ^ in the ID formalism, we pro- 
pose to conserve spin polarized current J-f/j,. Actually, 
we have used this boundary condition already in Sec- 
tion ^TOJ This boundary condition is justified when (i) 
thickness of the layer Li is thick enough compared to 
the change of pillar diameter between neighboring lay- 
ers, \a i+1 - cii\ < Li, \di - ai_ x \ < L { , where a, « y^, 
so that current has enough space to spread to different 
cross-sectional area Si (ii) pillar diameter is smaller 
than spin-flip length A. It can be shown that if Si — > inf , 
then j-m, /^T/i' MR etc. converge. So, it is not impor- 
tant, how large value Si is used to describe infinitely large 
layers. 

Fi sure [T^l shows j S p and fiy i calculated by 3D models 
in the case of constriction type structure. These data are 
compared with ID models with and without expanding 
the terminating Cu layers, called "ID constriction" and 
"S constant", respectively. 

3D constriction is described by ID constriction (S con- 
stant) model with precision of 15% (50%) for j sp and 
40% (80%) for An [Fig. H2» (b)] . The agreement be- 
tween Afj, for constriction type structure is quite poor, 
because Afi at this configuration is very small. In the 
case of Afj, larger than the above case, the agreement is 
about 10%. 

This ID model has been used to calculate j sp , Afi and 
MR in the position of free Co2 layer for all types of stud- 
ied structures. The results are presented in Table [I] in 
square brackets; we can see rather good agreement in 
all cases. An exception is the case of constriction type 
structure with infinitely large Col layer, as in this case 
the condition (ii) is not fulfilled. 



VI. CONCLUSION 

We have developed formalism which allows to calculate 
spin-polarized current j sp , and electrochemical potential 
fJtyi inside arbitrary electric circuit, consisting of ferro- 
or non-magnetic metallic SDRE elements as well as inter- 
face and surface resistivities. The formalism is limited to 
the parallel/antiparallel magnetic orientation in diffusive 
regime. 

To calculate spatial distribution of Hyi, jyi inside 
nanostructure, we divide the structure into an ID, 2D 
or 3D electric circuit network of SDRE elements which 
is successively solved. When division is carried out as 
described in Section IIII Al the renormalization of spin- 
flip length A has to be performed. 

This formalism is applied to (Cu/Co) 2 and (Cu/Co) 3 , 
pillar structures, where pillar cross-sectional area of 
starting/terminating layers were assumed to be either in- 
finitely large (column type, constriction type) or they 
have the same cross-sectional area as pillar (S const 
types). 

Inside the pillar surrounded by infinitely large layers, 
the j S p, j c h can be inhomogeneous. Maximal inhomo- 
geneity is found to be 29% and 20% in the case of in- 
finitely large Col layer. On the other hand, profile of 
Afj is found more homogeneous, with found maximum 
inhomogeneity of 2%. Such j sp inhomogeneities may lo- 
cally enhance the value of j sp , but they may disturb the 
magnetization reversal associated with spin precession. 
Due to oco < c C u, ARa.Co >> ARa,C U) both j ch and 
jsp flow rather perpendicular to the Co layers and fur- 
thermore the presence of Co layer makes j sp , j c h more 
homogeneous. 

When pillar is terminated by infinitely large layers, 
they serve as spin-scatterers, shortcutting up and down 
channels and hence modifying profiles of Afj and j S p- 

When such a spin-scatterers (but it is also valid for 
different spin-scatterers, as layers of Au, Ag, Pt, Ru, etc.) 
are introduced bellow "fixed" Col layer, they make Afj 
source "harder" . When they are placed above free Co2 
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layer, they shortcut up and down channels, reducing A/i 
nearly to zero, and hence enhancing j sp . 

Consequently, to get maximum j sp in the case of 
(Cu/Co) 2 at the position of "free" Co2 layer, it is impor- 
tant to introduce spin-scatterers both bellow fixed Col 
layer and above free Co2 layer. To get maximum A/i, it 
is important to introduce spin-scatterer bellow fixed Col 
layer, but reduce spin-scattering above Co2 free layer, 
latter one can be realized by reducing volume of material 
above Co2 layer, for example current drain can be thin 
long nanowire or cover layer with very thin thickness. 
As most of experimentally studied (Cu/Co) 2 structures 
have infinitely large layers at both ends, A/i inside them 
is nearly zero. 



For (Cu/Co) 3 system, maximum j sp (A/x) is for par- 
allel (antiparallel) magnetization of first, last Co layer. 
When applying above described optimizations, j sp and 
A/i can be further enhanced. 

Furthermore, in agreement with experimental results 
30, 3lL l32j| we found that j sp is linearly proportional to 
AB, = Ri — i£j (and hence MR is increased when increas- 
ing jsp at the position of free Co2 layer). Furthermore, 
dependence of Afi(AR) is also linear, but A/x is reduced 
when increasing AR (when increasing MR). 

Finally, we propose simple modification of ID Valet- 
Fert formalism [ll|. to incorporate spin-scattering in- 
duced by infinitely large layers attached to pillar. 
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